Bilateral stochastic power grid dispatching method

ABSTRACT

The disclosure relates to a two-side stochastic dispatching method for a power grid. By analyzing historical data of wind power, the Gaussian mixture distribution is fitted by software. For certain power system parameters, a two-side chance-constrained stochastic dispatching model is established. The hyperbolic tangent function is used to analyze and approximate cumulative distribution functions of random variables in the reserve demand constraint and the power flow constraint, to convert the two-side chance constraint into a deterministic constraint. The disclosure can have the advantage of using the hyperbolic tangent function to convert the two-side chance constraint containing risk levels and random variables into the solvable deterministic convex constraint, effectively improving the solution efficiency of the model, and providing decision makers with a more reasonable dispatching basis.

CROSS-REFERENCE TO RELATED APPLICATION

This application claims priority to Chinese Patent Application No. 202110301505.1, filed Mar. 22, 2021, the entire disclosure of which is incorporated herein by reference.

TECHNICAL FIELD

The disclosure relates to the field of operation technologies for power systems, and more particularly to a two-side stochastic dispatching method for a power grid.

BACKGROUND

It is major measures for energy development strategy to develop and utilize renewable energy and to realize sustainable energy development. With the large-scale integration of renewable energy such as wind power and photovoltaic power into the power grid, volatility and randomness of the renewable energy have brought two problems to active power rolling dispatch of the power system.

On the one hand, accurate and flexible wind power forecasting is the basis for realizing safe and economical active power rolling dispatch. The conventional forecasting methods include an interval description method of given upper and lower power limits and a description method of simple Gaussian probability density function. Although models such as Beta distribution and general distribution are also used in the fitting of the predicted power of renewable energy, these models either may not accurately fit the predicted power of wind power, or bring great difficulties to the solution of the active power rolling dispatch model. Therefore, an accurate and flexible forecasting model may be in urgent need of application.

On the other hand, the volatility and randomness of renewable energy make it difficult to apply the conventional deterministic or robust dispatching methods. Single-side chance-constrained dispatching may take into account the system operation risk and the power generation cost, but the independent modeling of the upper and lower bound constrain of the reserve capacity (branch capacity) results in the optimization result being too slack and unable to meet the preset safety level. Two-side chance-constrained stochastic dispatching may overcome the disadvantages of the unilateral model and obtain the lowest-cost dispatching strategy by minimizing the objective function value. The solution of two-side chance-constrained stochastic dispatching is extremely difficult. The conventional sampling solution method has a large amount of calculation and poor convergence. The slack solution method makes the solution result not accurate enough to achieve the high efficiency of dispatching.

In summary, the modeling and rapid solution of two-side chance-constrained dispatching taking into account the randomness of renewable energy power is still a major problem affecting the utilization of renewable energy.

SUMMARY

The disclosure aims to provide a two-side stochastic dispatching method for a power grid. A two-side chance constraint is transformed into a deterministic convex constraint and the deterministic convex constraint is solved effectively based on a Gaussian mixture distribution being accurately approximated by a hyperbolic tangent function, thus making full use of advantages of chance-constrained stochastic dispatching, effectively reducing system risks and saving power grid dispatching cost.

The two-side stochastic dispatching method for the power grid, provided in the disclosure, includes the following steps:

(1) establishing a two-side stochastic dispatching model based on hyperbolic tangent function approximation, in which the two-side stochastic dispatching model includes an objective function and constraints;

(2) analyzing and approximating a Gaussian mixture distribution satisfied by a random variable, appearing in the constraints, based on a reserve demand constraint and a power flow constraint in step (1), to obtain an approximate expression of the random variable; (3) transforming the two-side stochastic dispatching model in step (1) into a deterministic optimization model based on the approximation of the random variable in the reserve demand constraint and the power flow constraint in step (2); and

(4) using an interior point method to solve the deterministic optimization model transformed in steps (1) to (3), to obtain a power generation plan of thermal power units and renewable energy power stations.

The characteristics and advantages of the two-side stochastic dispatching method based on hyperbolic tangent function approximation, provided in the disclosure, may be as follows.

The method of the disclosure first accurately depicts characteristics and correlation of the predicted power of renewable energy such as wind power/photovoltaic power through the Gaussian mixture distribution of multiple random variables. Based on the distribution, the method establishes the cost-minimized stochastic dispatching model considering deterministic constraints and two-side chance constraints. The two-side chance constraints limit the uncertainty of the reserve demand (branch flow), caused by the randomness of the power of renewable energy power stations such as wind power/photovoltaic power stations during the dispatching process, to be within the upper limit and lower limit of safety by a certain level of confidence. At the same time, the hyperbolic tangent function is used to analyze and approximate the cumulative distribution function of the random variable in the chance constraint, so as to convert the two-side chance constraint into the deterministic convex constraint. The stochastic dispatching model is effectively solved. The result of model optimization is the optimal dispatching decision of conventional thermal power units and renewable energy power stations such as wind power/photovoltaic power under the control of dispatching risk and the reduction of dispatching cost. The advantage of the disclosure is transforming the two-side chance constraint into the deterministic convex constraint through the analytical approximation of the cumulative distribution function of the random variable, which effectively improves the solution efficiency of the model. Also, the two-side chance-constrained model with the adjustable risk level eliminates the conservative nature of conventional robust dispatching, which provides decision-makers with a more reasonable dispatching basis. The method of the disclosure may be applied to the dispatching of the power system including large-scale renewable energy grid connection.

DETAILED DESCRIPTION

The two-side stochastic dispatching method based on hyperbolic tangent function approximation, provided in the disclosure, includes the following steps:

The two-side stochastic dispatching method for the power grid, provided in the disclosure, includes the following steps:

(1) establishing a two-side stochastic dispatching model based on hyperbolic tangent function approximation, in which the two-side stochastic dispatching model includes an objective function and constraints;

(2) analyzing and approximating a Gaussian mixture distribution satisfied by a random variable, appearing in the constraints, based on a reserve demand constraint and a power flow constraint in step (1), to obtain an approximate expression of the random variable; (3) transforming the two-side stochastic dispatching model in step (1) into a deterministic optimization model based on the approximation of the random variable in the reserve demand constraint and the power flow constraint in step (2); and

(4) using an interior point method to solve the deterministic optimization model transformed in steps (1) to (3), to obtain a power generation plan of thermal power units and renewable energy power stations.

In the above two-side stochastic dispatching method for the power grid, step (1) includes:

(1) establishing the objective function of the two-side stochastic dispatching model, including:

minimizing a sum of power generation cost, positive spinning reserve capacity cost and negative spinning reserve capacity cost, of the thermal power units, which is expressed by a formula of:

$\begin{matrix} {{\min{\sum\limits_{t = 1}^{T}{\sum\limits_{i = 1}^{N_{G}}{{CF}_{i}\left( P_{i}^{t} \right)}}}} + {\sum\limits_{t = 1}^{T}{\sum\limits_{i = 1}^{N_{G}}\left( {{c_{up}r_{i,t}^{up}} + {c_{dn}r_{i,t}^{dn}}} \right)}}} & (25) \end{matrix}$

where, T represents a number of dispatching periods,

N_(G) represents a number of the thermal power units,

t represents a serial number of the dispatching period,

i represents a serial number of the thermal power unit,

c_(up) represents a positive spinning reserve capacity cost per unit,

c_(dn) represents a negative spinning reserve capacity cost per unit,

r_(i,t) ^(up) represents a positive spinning reserve capacity reserved by an i^(th) thermal power unit during a dispatching period t,

r_(i,t) ^(dn) represents a negative spinning reserve capacity reserved by the i^(th) thermal power unit during the dispatching period t,

P_(i) ^(t) represents a planned power of the i^(th) thermal power unit during the dispatching period t,

CF_(i) represents a fuel cost function of the i^(th) thermal power unit,

the fuel cost function of the thermal power unit is expressed as a quadratic function of the power P_(i) ^(t) of the thermal power unit by a formula of: CF_(i)(P _(i) ^(t))=a _(i)(P _(i) ^(t))² +b _(i) P _(i) ^(t) +c _(i)  (26)

where, a_(i) represents a quadratic coefficient of a fuel cost of the i^(th) thermal power unit,

b_(i) represents a primary coefficient of the fuel cost of the i^(th) thermal power unit,

c_(i) represents a constant coefficient of the fuel cost of the i^(th) thermal power unit,

(2) the constraints of the two-side stochastic dispatching model, include:

(2-1) a power balance constraint of the power grid, which is expressed by a formula of:

$\begin{matrix} {{{\sum\limits_{i = 1}^{N_{G}}P_{i}^{t}} + {\overset{N_{W}}{\sum\limits_{j = 1}}w_{j}^{t}}} = {\sum\limits_{k = 1}^{N_{D}}d_{k}^{t}}} & (27) \end{matrix}$

where, w^(t) _(j) represents a planned reference power of a j^(th) renewable energy power station during the dispatching period t,

d^(t) _(k) represents a size of a k^(th) load during the dispatching period t,

N_(w) represents a number of the power stations, and

N_(D) represents a number of loads;

(2-2) an upper limit constraint and a lower limit constraint of the power P_(i) ^(t) of the thermal power unit in the power grid, which is expressed by a formula of: P _(i) ≤P _(i) ^(t) ≤P _(i)  (28)

where, P _(i) represents an upper bound of the power of the i^(th) thermal power unit, and

P _(i) represents a lower bound of the power of the i^(th) thermal power unit;

(2-3) an upper limit constraint and a lower limit constraint of the planned power w^(t) _(j) of the power station of the power grid, which is expressed by a formula of: 0≤w ^(t) _(j) ≤w ^(t) _(j)  (29)

where, w ^(t) _(j) represents an upper bound of an allowable power of the j^(th) renewable energy power station during the dispatching period t;

(2-4) ramp constraints of the thermal power unit in the power grid, which are expressed by formulas of: P _(i) ^(t) −P _(i) ^(t−1)≥−RD_(i) ΔT  (30) and P _(i) ^(t) −P _(i) ^(t−1)≤RU_(i) ΔT  (31)

where, RU_(i) represents an upward ramp rate of the i^(th) thermal power unit,

RD_(i) represents a downward ramp rate of the i^(th) thermal power unit,

ΔT represents a dispatching interval between two adjacent dispatching periods;

(2-5) a maximum reserve capacity constraint that the thermal power unit in the power grid provides, which is expressed by formulas of: r _(i,t) ^(up)≤min{ P _(i) −P _(i) ^(t),RU_(i) ΔT}  (32) and r _(i,t) ^(dn)≤min{P _(i) ^(t) −P _(i),RD_(i) ΔT}  (33)

where, a mathematical symbol min {·} represents a minimum value of elements in a set;

(2-6) a reserve demand constraint in response to fluctuations in the renewable energy power stations in the power grid, which is expressed by a formula of:

Pr ( - ∑ i = 1 N G r i , t up ≤ ∑ j = 1 N W w ~ j t - ∑ j = 1 N W w j t ≤   ∑ i = 1 N G r i , t d ⁢ n ) ≥ 1 - β ( 34 )

where, a mathematical symbol Pr{·} represents a probability of occurrence of an event,

{tilde over (w)}^(t) _(j) represents an actual power of the j^(th) renewable energy power station during the dispatching period t,

a superscript ˜ represents that a variable is a random variable,

β represents a maximum risk level preset by a dispatcher and a value range of the maximum risk level is 0 to 0.5, a joint probability distribution of the actual power of all renewable energy power stations meets a Gaussian mixture distribution of:

$\begin{matrix} {{{\overset{\sim}{w}}^{t} = {\left\{ {\overset{\sim}{w}}_{j}^{t} \middle| {1 \leq j \leq N_{W}} \right\}{and}}}{{{{PDF}_{{\overset{\sim}{w}}^{t}}(Y)} = {\sum\limits_{m = 1}^{M}{\omega_{m}{N\left( {Y,\mu_{m},\Sigma_{m}} \right)}}}},{\omega_{m} \geq 0}}} & (35) \end{matrix}$ $\begin{matrix} {{N\left( {\left. Y \middle| \mu_{m} \right.,\Sigma_{m}} \right)} = {\frac{1}{\left( {2\pi} \right)^{N_{W}/2}{\det\left( \Sigma_{m} \right)}^{1/2}}e^{{- \frac{1}{2}}{({Y - \mu_{m}})}^{T}{\Sigma_{m}^{- 1}({Y - \mu_{m}})}}}} & (36) \end{matrix}$

{tilde over (w)}^(t) represents a planned power set of all renewable energy power stations during the dispatching period t,

{tilde over (w)}^(t) is a stochastic vector,

PDF_({tilde over (w)}) _(t) (·) represents a probability density function of the stochastic vector,

Y represents a value of {tilde over (w)}^(t),

N(Y, μ_(m), Σ_(m)) represents a m^(th) component of the Gaussian mixture distribution,

M represents a number of components of the Gaussian mixture distribution,

ω_(m) represents a weight coefficient of the m^(th) component of the Gaussian mixture distribution and satisfies that a sum of weight coefficients of all components is equal to 1,

μ_(m) represents a mean vector of the m^(th) component,

Σ_(m) represents a covariance matrix of the m^(th) component,

det represents a determinant of the matrix, and

a superscript T represents a transpose of the matrix;

(2-7) a power flow constraint of the power grid, which is expressed by a formula of:

$\begin{matrix} {{{P{r\left( {{- L_{l}} \leq {{\sum\limits_{i = 1}^{N_{G}}{G_{l,i}{\overset{\sim}{P}}_{i}^{t}}} + {\sum\limits_{j = 1}^{N_{W}}{G_{l,j}{\overset{\sim}{w}}_{j}^{t}}} - {\sum\limits_{k = 1}^{N_{D}}{G_{l,k}d_{k}^{t}}}} \leq L_{l}} \right)}} \geq {1 - \eta}},} & (37) \end{matrix}$

where,

G_(l,i) represents a transfer distribution factor of a l^(th) branch for an active power of the i^(th) thermal power unit,

G_(l,j) represents a transfer distribution factor of a l^(th) branch for an active power of the j^(th) renewable energy power station,

G_(l,k) represents a transfer distribution factor of the l^(th) branch for the k^(th) load,

the above-mentioned transfer distribution factors are respectively obtained from a power grid dispatching center,

L_(l) represents an upper limit of active power on the l^(th) branch,

η represents a risk level that the active power on the branch exceeds the upper limit of active power on the branch, which is set by the dispatcher,

{tilde over (P)}_(i) ^(t) represents an actual power of the i^(th) thermal power unit during the dispatching period t, which is expressed by a formula of:

$\begin{matrix} {{{\overset{\sim}{P}}_{i}^{t} = {P_{i}^{t} - {\alpha_{i}\left( {{\sum\limits_{j = 1}^{N_{W}}{\overset{\sim}{w}}_{j}^{t}} - {\sum\limits_{j = 1}^{N_{W}}w_{j}^{t}}} \right)}}}{{\sum\limits_{i = 1}^{N_{G}}\alpha_{i}}\  = 1}} & (38) \end{matrix}$

where, α_(i) represents a power distribution coefficient of the i^(th) thermal power unit, which is determined by a ratio of a rated capacity of the thermal power unit to a total capacity of all thermal power units participating in power generation;

in order to meet the power balance of the power grid at any dispatch time, a sum of power distribution coefficients of all thermal Dower units is equal to 1,

${\sum\limits_{j = 1}^{N_{W}}{\overset{\sim}{w}}_{j}^{t}} - {\sum\limits_{j = 1}^{N_{W}}w_{j}^{t}}$ represents a power difference between actual powers and planned powers of all renewable energy power stations;

substituting the formula (14) into the formula (13), in which the formula (13) is transformed into a formula of:

$\begin{matrix} {{\Pr\left( {{- L_{l}}\  \leq {{\sum\limits_{i = 1}^{N_{G}}{G_{l,i}P_{i}^{t}}} + {\sum\limits_{i = 1}^{N_{G}}{G_{l,i}\alpha_{i}{\sum\limits_{j = 1}^{N_{W}}w_{j}^{t}}}} + {\sum\limits_{j = 1}^{N_{W}}{\left( {G_{l,j} - {\sum\limits_{i = 1}^{N_{G}}{G_{l,i}\alpha_{i}}}} \right){\overset{\sim}{w}}_{j}^{t}}} - {\sum\limits_{k = 1}^{N_{D}}{G_{l,k}d_{k}^{t}}}} \leq L_{l}} \right)} \geq {1 - {\eta.}}} & (39) \end{matrix}$

In the above two-side stochastic dispatching method for the power grid, step (2) includes:

(1) approximating the Gaussian mixture distribution by a formula of:

$\begin{matrix} {{{CD{F_{\xi_{0}}(x)}} \approx {\sum\limits_{m_{0}}^{M_{0}}{\frac{\omega_{0}}{2}\left( {1 + {\tanh\left( {{a\frac{x - \mu_{m_{0}}}{\sigma_{m_{0}}}} + {b\left( \frac{x - \mu_{m_{0}}}{\sigma_{m_{0}}} \right)}^{3}} \right)}} \right)}}}{{a = {{0.7}983}},{b = {{0.0}3564}}}{{\sum\limits_{m_{0}}^{M_{0}}\omega_{0}} = 1}} & (40) \end{matrix}$

where, ξ₀ represents any one-dimensional random variable that obeys the Gaussian mixture distribution,

CDF_(ξ) ₀ represents a cumulative distribution function of the random variable,

x represents an independent variable of the cumulative distribution function and represents a possible value of ξ₀,

M₀ represents a number of Gaussian components obeyed by the random variable ξ₀,

m₀ represents a m₀ ^(th) component,

ω₀, μ_(m) ₀ and σ_(m) ₀ represent a weight coefficient, an expectation and a standard deviation of the m₀ ^(th) component,

a and b represent constants, and

tan h represent a hyperbolic tangent function in mathematical operations;

(2) based on the approximating method of the Gaussian mixture distribution in step (1), approximating a cumulative distribution function of a random variable in the reserve demand constraint in the formula (13) by a formula of:

$\begin{matrix} {{{{CDF}_{\xi_{r}^{t}}(y)} \approx {{ACDF}_{\xi_{r}^{t}}(y)}} = {\sum\limits_{m}^{M}{\frac{\omega_{m}}{2}\left( {1 + {\tanh\left( {{a\frac{y - {1^{T}\mu_{m}}}{1^{T}\Sigma_{m}1}} + {b\left( \frac{y - {1^{T}\mu_{m}}}{1^{T}\Sigma_{m}1} \right)}^{3}} \right)}} \right)}}} & (41) \end{matrix}$ ${\xi_{r}^{t} = {1^{T}{\overset{\sim}{w}}_{j}^{t}}},$

where, ξ_(r) ^(t) represents a random variable and represents a sum of actual powers of all renewable energy power stations at time t,

1 represents a column vector with all elements being 1,

CDF_(ξ) _(r) _(t) (·) represents a probability density function of the one-dimensional random variable ξ_(r) ^(t),

y represents an independent variable of the cumulative distribution function CDF_(ξ) _(r) _(t) (·) and represents a possible value of ξ_(r) ^(t), and

ACDF_(ξ) _(r) _(t) (·) represents an approximate expression of CDF_(ξ) _(r) _(t) (·); and

(3) based on the approximating method of the Gaussian mixture distribution in step (1), approximating a cumulative distribution function of a random variable in the power flow constraint of the power grid in the formula (15) by a formula of:

$\begin{matrix} {{{{CDF}_{\xi_{l}^{t}}(z)} \approx {{ACDF}_{\xi_{l}^{t}}(z)}} = {\sum\limits_{m}^{M}{\frac{\omega_{m}}{2}\left( {1 + {\tanh\left( {{a\frac{z - {F_{l}^{T}\mu_{m}}}{F_{l}^{T}\Sigma_{m}F_{l}}} + {b\left( \frac{z - {F_{l}^{T}\mu_{m}}}{F_{l}^{T}\Sigma_{m}F_{l}} \right)}^{3}} \right)}} \right)}}} & (42) \end{matrix}$ $\xi_{l}^{t} = {F_{l}^{T}{\overset{\sim}{w}}_{j}^{t}}$

where, F_(l) represents a column vector and a j^(th) element of the vector is

${G_{l,j} - {\sum\limits_{i = 1}^{N_{G}}{G_{l,i}\alpha_{i}}}},$

ξ_(l) ^(t) represents a random variable and represents a linear combination of actual powers of all renewable energy power stations at time t, coefficients of the linear combination is F_(l),

CDF_(ξ) _(l) _(t) (·) represents a probability density function of the one-dimensional random variable ξ_(l) ^(t),

z represents an independent variable of the cumulative distribution function CDF_(ξ) _(l) _(t) (·) and represents a possible value of and ξ_(l) ^(t), and

ACDF_(ξ) _(l) _(t) (·) represents an approximate expression of CDF_(ξ) _(l) _(t) (·).

In the above two-side stochastic dispatching method for the power grid, step (4) includes:

(1) transforming a chance constraint by using the approximate expression (17) of the cumulative distribution function based on the approximation of the random variable in the reserve demand constraint, as formulas of:

$\begin{matrix} {{{\sum\limits_{i = 1}^{N_{G}}r_{i,t}^{dn}} + {\sum\limits_{j = 1}^{N_{W}}w_{j}^{t}}} \geq q_{r,t}^{up}} & (43) \end{matrix}$ $\begin{matrix} {{{- {\sum\limits_{i = 1}^{N_{G}}r_{i,t}^{up}}} + {\sum\limits_{j = 1}^{N_{W}}w_{j}^{t}}} \leq q_{r,t}^{dn}} & (44) \end{matrix}$ $\begin{matrix} {{{{ACDF}_{\xi_{r}^{t}}\left( q_{r,t}^{up} \right)} - {{ACDF}_{\xi_{r}^{t}}\left( q_{r,t}^{dn} \right)}} \geq {1 - \beta}} & (45) \end{matrix}$

where, q_(r,t) ^(up) and q_(r,t) ^(dn) are auxiliary variables, which respectively represent a upper and lower quantiles of the reserve demand constraint at time t;

(2) transforming a power flow chance constraint determined in the formula (15) by using the approximate expression (18) of the cumulative distribution function of the random variable in the power flow constraint of the power grid based on the approximation of the random variable in the power flow constraint of the power grid, as formulas of:

$\begin{matrix} {{L_{l} - {\sum\limits_{i = 1}^{N_{G}}{G_{l,i}P_{i}^{t}}} - {\sum\limits_{i = 1}^{N_{G}}{G_{l,i}\alpha_{i}{\sum\limits_{j = 1}^{N_{W}}w_{j}^{t}}}} + {\sum\limits_{k = 1}^{N_{D}}{G_{l,k}d_{k}^{t}}}} \geq q_{l,t}^{up}} & (46) \end{matrix}$ $\begin{matrix} {{{- L_{l}} - {\sum\limits_{i = 1}^{N_{G}}{G_{l,i}P_{i}^{t}}} - {\sum\limits_{i = 1}^{N_{G}}{G_{l,i}\alpha_{i}{\sum\limits_{j = 1}^{N_{W}}w_{j}^{t}}}} + {\sum\limits_{k = 1}^{N_{D}}{G_{l,k}d_{k}^{t}}}} \leq q_{l,t}^{up}} & (47) \end{matrix}$ $\begin{matrix} {{{{ACDF}_{\xi_{l}^{t}}\left( q_{l,t}^{up} \right)} - {{ACDF}_{\xi_{l}^{t}}\left( q_{l,t}^{dn} \right)}} \geq {1 - \beta}} & (48) \end{matrix}$

where, q_(l,t) ^(up) and q_(l,t) ^(dn) are auxiliary variables, which respectively represent a upper and lower quantiles of the power flow constraint of the l^(th) branch at time t;

(4) solving the two-side stochastic dispatching model based on hyperbolic tangent function approximation, in which,

the interior point method is used to solve the optimization model including the objective function (1)-(2) and constraints (3)-(9), (14), (19)-(24) in the above steps (1)-(3), to obtain the planned power P_(i) ^(t) of the i^(th) thermal power unit during the dispatching period t and the planned reference power w^(t) _(j) of the j^(th) renewable energy power station during the dispatching period t, in which P_(i) ^(t) is regarded as the planned power of the i^(th) thermal power unit during the dispatching period t and w^(t) _(j) is regarded as the reference power of the j^(th) renewable energy power station during the dispatching period t, to realize the chance-constrained stochastic economic dispatching based on the approximation of the hyperbolic tangent function. 

What is claimed is:
 1. A two-side stochastic dispatching method for a power grid, comprising: (1) establishing, by the power grid, a two-side stochastic dispatching model based on a hyperbolic tangent function approximation, wherein the two-side stochastic dispatching model comprises an objective function and constraints; (2) analyzing and approximating, by the power grid, a Gaussian mixture distribution of a random variable appearing in the constraints, based on a reserve demand constraint and a power flow constraint of the constraints in step (1), to obtain an approximate expression of the random variable; (3) transforming, by the power grid, the two-side stochastic dispatching model in step (1) into a deterministic optimization model based on the approximation of the random variable in the reserve demand constraint and the power flow constraint in step (2); (4) using, by the power grid, an interior point method to solve the deterministic optimization model transformed in steps (1) to (3), to obtain a power generation plan of thermal power units and renewable energy power stations of the power grid; and operating, by the power grid, the thermal power units and the renewable energy power stations based on the power generation plan.
 2. The method according to claim 1, wherein the step (1) comprises: (1) the establishing the objective function of the two-side stochastic dispatching model, comprising: minimizing a sum of a power generation cost, a positive spinning reserve capacity cost and a negative spinning reserve capacity cost, of the thermal power units, which is expressed by a formula of: $\begin{matrix} {{\min{\sum\limits_{t = 1}^{T}{\sum\limits_{i = 1}^{N_{G}}{{CF}_{i}\left( P_{i}^{t} \right)}}}} + {\sum\limits_{t = 1}^{T}{\sum\limits_{i = 1}^{N_{G}}\left( {{c_{up}r_{i,t}^{up}} + {c_{dn}r_{i,t}^{dn}}} \right)}}} & (1) \end{matrix}$ where, T represents a number of dispatching periods, N_(G) represents a number of the thermal power units, t represents a serial number of the dispatching period, i represents a serial number of the thermal power unit, c_(up) represents a positive spinning reserve capacity cost per unit, c_(dn) represents a negative spinning reserve capacity cost per unit, r_(i,t) ^(up) represents a positive spinning reserve capacity reserved by an i^(th) thermal power unit during a dispatching period t, r_(i,t) ^(dn) represents a negative spinning reserve capacity reserved by the i^(th) thermal power unit during the dispatching period t, P_(i) ^(t) represents a planned power of the i^(th) thermal power unit during the dispatching period t, CF_(i) represents a fuel cost function of the i^(th) thermal power unit, wherein the fuel cost function of the thermal power unit is expressed as a quadratic function of the power P_(i) ^(t) of the thermal power unit by a formula of: CF_(i)(P _(i) ^(t))=a _(i)(P _(i) ^(t))² +b _(i) P _(i) ^(t) +c _(i)  (2) where, a_(i) represents a quadratic coefficient of a fuel cost of the i^(th) thermal power unit, b_(i) represents a primary coefficient of the fuel cost of the i^(th) thermal power unit, c_(i) represents a constant coefficient of the fuel cost of the i^(th) thermal power unit, (2) the constraints of the two-side stochastic dispatching model, comprising: (2-1) a power balance constraint of the power grid, which is expressed by a formula of: $\begin{matrix} {{{\sum\limits_{i = 1}^{N_{G}}P_{i}^{t}} + {\overset{N_{W}}{\sum\limits_{j = 1}}w_{j}^{t}}} = {\sum\limits_{k = 1}^{N_{D}}d_{k}^{t}}} & (3) \end{matrix}$ where, w_(j) ^(t) represents a planned reference power of a j^(th) renewable energy power station during the dispatching period t, d_(k) ^(t) represents a size of a k^(th) load during the dispatching period t, N_(w) represents a number of the power stations, and N_(D) represents a number of loads; (2-2) an upper limit constraint and a lower limit constraint of the power P_(i) ^(t) of the thermal power unit in the power grid, which is expressed by a formula of: P _(i) ≤P _(i) ^(t) ≤P _(i)  (4) where, P _(i) represents an upper bound of the power of the i^(th) thermal power unit, and P _(i) represents a lower bound of the power of the i^(th) thermal power unit; (2-3) an upper limit constraint and a lower limit constraint of the planned power w^(t) _(j) of the power station of the power grid, which is expressed by a formula of: 0≤w ^(t) _(j) ≤w ^(t) _(j)  (5) where, w ^(t) _(j) represents an upper bound of an allowable power of the j^(th) renewable energy power station during the dispatching period t; (2-4) ramp constraints of the thermal power unit in the power grid, which are expressed by formulas of: P _(i) ^(t) −P _(i) ^(t−1)≥−RD_(i) ΔT  (6) and P _(i) ^(t) −P _(i) ^(t−1)≤RU_(i) ΔT  (31) where, RU_(i) represents an upward ramp rate of the i^(th) thermal power unit, RD_(i) represents a downward ramp rate of the i^(th) thermal power unit, ΔT represents a dispatching interval between two adjacent dispatching periods; (2-5) a maximum reserve capacity constraint that the thermal power unit in the power grid provides, which is expressed by formulas of: r _(i,t) ^(up)≤min{ P _(i) −P _(i) ^(t),RU_(i) ΔT}  (8) and r _(i,t) ^(dn)≤min{P _(i) ^(t) −P _(i),RD_(i) ΔT}  (9) where, a mathematical symbol min {·} represents a minimum value of elements in a set; (2-6) the reserve demand constraint in response to fluctuations in the renewable energy power stations in the power grid, which is expressed by a formula of: $\begin{matrix} {{P{r\left( {{- {\sum\limits_{i = 1}^{N_{G}}r_{i,t}^{up}}} \leq {{\overset{N_{W}}{\sum\limits_{j = 1}}{\overset{\sim}{w}}_{j}^{t}} - {\overset{N_{W}}{\sum\limits_{j = 1}}w_{j}^{t}}} \leq {\sum\limits_{i = 1}^{N_{G}}r_{i,t}^{dn}}} \right)}} \geq {1 - \beta}} & (10) \end{matrix}$ where, a mathematical symbol Pr {·} represents a probability of occurrence of an event, {tilde over (w)}^(t) _(j) represents an actual power of the j^(th) renewable energy power station during the dispatching period t, a superscript ˜ represents that a variable is a random variable, β represents a maximum risk level preset by a dispatcher and a value range of the maximum risk level is 0 to 0.5, a joint probability distribution of the actual power of all of the renewable energy power stations meets a Gaussian mixture distribution of: $\begin{matrix} {{\overset{\sim}{w}}^{t} = {\left\{ {\overset{\sim}{w}}_{j}^{t} \middle| {1 \leq j \leq N_{W}} \right\}{and}}} & (11) \end{matrix}$ $\begin{matrix} {{{{{PDF}_{{\overset{\sim}{w}}^{t}}\ (Y)} = {\sum\limits_{m = 1}^{M}{\omega_{m}{N\left( {Y,\mu_{m},\sum_{m}} \right)}}}},{\omega_{m} \geq 0}}{{N\left( {\left. Y \middle| \mu_{m} \right.,\sum_{m}} \right)} = {\frac{1}{\left( {2\pi} \right)^{N_{W}/2}{\det\left( \sum_{m} \right)}^{1/2}}e^{{- \frac{1}{2}}{({Y - \mu_{m}})}^{T}{\sum_{m}^{- 1}{({Y - \mu_{m}})}}}}}} & (12) \end{matrix}$ {tilde over (w)}^(t) represents a planned power set of all of the renewable energy power stations during the dispatching period t, {tilde over (w)}^(t) is a stochastic vector, PDF_({tilde over (w)}) _(t) (·) represents a probability density function of the stochastic vector, Y represents a value of {tilde over (w)}^(t), N(Y, μ_(m), Σ_(m)) represents a m^(th) component of the Gaussian mixture distribution, M represents a number of components of the Gaussian mixture distribution, ω_(m), represents a weight coefficient of the m^(th) component of the Gaussian mixture distribution and satisfies that a sum of weight coefficients of all components is equal to 1, μ_(m), represents a mean vector of the m^(th) component, Σ_(m) represents a covariance matrix of the m^(th) component, det represents a determinant of the matrix, and a superscript T represents a transpose of the matrix; (2-7) the power flow constraint of the power grid, which is expressed by a formula of: $\begin{matrix} {{{P{r\left( {{- L_{l}} \leq {{\sum\limits_{i = 1}^{N_{G}}{G_{l}{\overset{\sim}{P}}_{i}^{t}}} + {\overset{N_{W}}{\sum\limits_{j = 1}}{G_{l,j}{\overset{\sim}{w}}_{j}^{t}}} - {\sum\limits_{k = 1}^{N_{D}}{G_{l,k}d_{k}^{t}}}} \leq L_{l}} \right)}} \geq {1 - \eta}},} & (13) \end{matrix}$ where, G_(l,j) represents a transfer distribution factor of a l^(th) branch for an active power of the i^(th) thermal power unit, G_(l,j) represents a transfer distribution factor of the l^(h) branch for an active power of the j^(th) renewable energy power station, G_(l,k) represents a transfer distribution factor of the l^(th) branch for the k^(th) load, the above-mentioned transfer distribution factors are respectively obtained from a power grid dispatching center, L_(l) represents an upper limit of active power on the l^(th) branch, η represents a risk level that the active power on the branch exceeds the upper limit of active power on the branch, which is set by the dispatcher, {tilde over (P)}_(i) ^(t) represents an actual power of the i^(th) thermal power unit during the dispatching period t, which is expressed by a formula of: $\begin{matrix} {{{\overset{\sim}{P}}_{i}^{t} = {P_{i}^{t} - {\alpha_{i}\left( {{\overset{N_{W}}{\sum\limits_{j = 1}}{\overset{\sim}{w}}_{j}^{t}} - {\overset{N_{W}}{\sum\limits_{j = 1}}w_{j}^{t}}} \right)}}}{{\sum\limits_{i = 1}^{N_{G}}\alpha_{i}} = 1}} & (14) \end{matrix}$ where, α_(i) represents a power distribution coefficient of the i^(th) thermal power unit, which is determined by a ratio of a rated capacity of the thermal power unit to a total capacity of all thermal power units participating in power generation; in order to meet the power balance of the power grid at any dispatch time, a sum of power distribution coefficients of all thermal power units is equal to 1, ${\overset{N_{W}}{\sum\limits_{j = 1}}{\overset{\sim}{w}}_{j}^{t}} - {\overset{N_{W}}{\sum\limits_{j = 1}}w_{j}^{t}}$  represents a power difference between actual powers and planned powers of all of the renewable energy power stations; substituting the formula (14) into the formula (13), wherein the formula (13) is transformed into a formula of: $\begin{matrix} {{P{r\left( {{- L_{l}} \leq {{\sum\limits_{i = 1}^{N_{G}}{G_{l,i}P_{i}^{t}}} + {\sum\limits_{i = 1}^{N_{G}}{G_{l,i}\alpha_{i}{\overset{N_{W}}{\sum\limits_{j = 1}}w_{j}^{t}}}} + {\overset{N_{W}}{\sum\limits_{j = 1}}{\left( {G_{l,j} - {\sum\limits_{i = 1}^{N_{G}}{G_{l,i}\alpha_{i}}}} \right){\overset{\sim}{w}}_{j}^{t}}} - {\sum\limits_{k = 1}^{N_{D}}{G_{l,k}d_{k}^{t}}}} \leq L_{l}} \right)}} \geq {{1 - \eta}.}} & (15) \end{matrix}$
 3. The method according to claim 1, wherein the step (2) comprises: (1) approximating the Gaussian mixture distribution by a formula of: $\begin{matrix} {{{CD{F_{\xi_{0}}(x)}} \approx {\sum\limits_{m_{0}}^{M_{0}}{\frac{\omega_{0}}{2}\left( {1 + {\tanh\left( {{a\frac{x - \mu_{m_{0}}}{\sigma_{m_{0}}}} + {b\left( \frac{x - \mu_{m_{0}}}{\sigma_{m_{0}}} \right)}^{3}} \right)}} \right)}}}{{a = {{0.7}983}},{b = {{0.0}3564}}}{{\sum\limits_{m_{0}}^{M_{0}}\omega_{0}} = 1}} & (16) \end{matrix}$ where, ξ₀ represents any one-dimensional random variable that obeys the Gaussian mixture distribution, CDF_(ξ) ₀ represents a cumulative distribution function of the random variable, x represents an independent variable of the cumulative distribution function and represents a possible value of ξ₀, M₀ represents a number of Gaussian components obeyed by the random variable ξ₀, m₀ represents a m₀ ^(th) component, ω₀, μ_(m) ₀ , σ_(m) ₀ represent a weight coefficient, an expectation and a standard deviation of the m₀ ^(th) component, a and b represent constants, and tanh represent a hyperbolic tangent function in mathematical operations; (2) based on the approximating method of the Gaussian mixture distribution in step (1), approximating the cumulative distribution function of the random variable in the reserve demand constraint in the formula (13) by a formula of: $\begin{matrix} {{{CD{F_{\xi_{r}^{t}}(y)}} \approx {ACD{F_{\xi_{r}^{t}}(y)}}} = {\sum\limits_{m}^{M}{\frac{\omega_{m}}{2}\left( {1 + {\tanh\left( {{a\frac{y - {1^{T}\mu_{m}}}{1^{T}{\sum_{m}1}}} + {b\left( \frac{y - {1^{T}\mu_{m}}}{1^{T}{\sum_{m}1}} \right)}^{3}} \right)}} \right)}}} & (17) \end{matrix}$ ${\xi_{r}^{t} = {1^{T}{\overset{\sim}{w}}_{j}^{t}}},$ where ξ_(r) ^(t) represents a random variable and represents a sum of actual powers of all of the renewable energy power stations at time t, 1 represents a column vector with all elements being 1, CDF_(ξ) _(r) _(t) (·) represents a probability density function of the one-dimensional random variable ξ_(r) ^(t), y represents an independent variable of the cumulative distribution function CDF_(ξ) _(r) _(t) (·) and represents a possible value of ξ_(r) ^(t), and ACDF_(ξ) _(r) _(t) (·) represents an approximate expression of CDF_(ξ) _(r) _(t) (·); and (3) based on the approximating method of the Gaussian mixture distribution in step (1), approximating the cumulative distribution function of the random variable in the power flow constraint of the power grid in the formula (15) by a formula of: $\begin{matrix} {{{CD{F_{\xi_{l}^{t}}(z)}} \approx {ACD{F_{\xi_{l}^{t}}(z)}}} = {\sum\limits_{m}^{M}{\frac{\omega_{m}}{2}\left( {1 + {\tanh\left( {{a\frac{z - {F_{l}^{T}\mu_{m}}}{F_{l}^{T}{\sum_{m}F_{l}}}} + {b\left( \frac{z - {F_{l}^{T}\mu_{m}}}{F_{l}^{T}{\sum_{m}F_{l}}} \right)}^{3}} \right)}} \right)}}} & (18) \end{matrix}$ $\xi_{l}^{t} = {F_{l}^{T}{\overset{\sim}{w}}_{j}^{t}}$ where, F_(l) represents a column vector and a j^(th) element of the vector is ${G_{l,j} - {\sum\limits_{i = 1}^{N_{G}}{G_{l,i}\alpha_{i}}}},$ ξ_(l) ^(t) represents a random variable and represents a linear combination of actual powers of all of the renewable energy power stations at time t, coefficients of the linear combination is F_(l), CDF_(ξ) _(l) _(t) (·) represents a probability density function of the one-dimensional random variable ξ_(l) ^(t), z represents an independent variable of the cumulative distribution function CDF_(ξ) _(l) _(t) (·) and represents a possible value of ξ_(l) ^(t), and ACDF_(ξ) _(l) _(t) (·) represents an approximate expression of CDF_(ξ) _(l) _(t) (·).
 4. The method according to claim 1, wherein the step (4) comprises: (1) transforming a chance constraint by using the formula (17) of the cumulative distribution function based on the approximation of the random variable in the reserve demand constraint, as formulas of: $\begin{matrix} {{{\sum\limits_{i = 1}^{N_{G}}r_{i,t}^{dn}} + {\sum\limits_{j = 1}^{N_{W}}w_{j}^{t}}} \geq q_{r,t}^{up}} & (19) \end{matrix}$ $\begin{matrix} {{{- {\sum\limits_{i = 1}^{N_{G}}r_{i,t}^{dn}}} + {\sum\limits_{j = 1}^{N_{W}}w_{j}^{t}}} \leq q_{r,t}^{up}} & (20) \end{matrix}$ $\begin{matrix} {{{{ACDF}_{\xi_{r}^{t}}\left( q_{r,t}^{up} \right)} - {ACD{F_{\xi_{r}^{t}}\left( q_{r,t}^{dn} \right)}}} \geq {1 - \beta}} & (21) \end{matrix}$ where, q_(r,t) ^(up) and q_(r,t) ^(dn) are auxiliary variables, which respectively represent a upper and lower quantiles of the reserve demand constraint at time t; (2) transforming a power flow chance constraint determined in the formula (15) by using the formula (18) of the cumulative distribution function of the random variable in the power flow constraint of the power grid based on the approximation of the random variable in the power flow constraint of the power grid, as formulas of: $\begin{matrix} {{L_{l} - {\sum\limits_{i = 1}^{N_{G}}{G_{l,i}P_{i}^{t}}} - {\sum\limits_{i = 1}^{N_{G}}{G_{l,i}\alpha_{i}{\sum\limits_{j = 1}^{N_{W}}w_{j}^{t}}}} + {\sum\limits_{k = 1}^{N_{D}}{G_{l,k}d_{k}^{t}}}} \geq q_{l,t}^{up}} & (22) \end{matrix}$ $\begin{matrix} {{{- L_{l}} - {\sum\limits_{i = 1}^{N_{G}}{G_{l,i}P_{i}^{t}}} - {\sum\limits_{i = 1}^{N_{G}}{G_{l,i}\alpha_{i}{\sum\limits_{j = 1}^{N_{W}}w_{j}^{t}}}} + {\sum\limits_{k = 1}^{N_{D}}{G_{l,k}d_{k}^{t}}}} \leq q_{l,t}^{up}} & (23) \end{matrix}$ $\begin{matrix} {{{{ACDF}_{\xi_{l}^{t}}\left( q_{l,t}^{up} \right)} - {{ACDF}_{\xi_{l}^{t}}\left( q_{l,t}^{dn} \right)}} \geq {1 - \beta}} & (24) \end{matrix}$ where, q_(l,t) ^(up) and q_(l,t) ^(dn) are auxiliary variables, which respectively represent a upper and lower quantiles of the power flow constraint of the l^(th) branch at time t; (4) solving the two-side stochastic dispatching model based on the hyperbolic tangent function approximation, wherein, the interior point method is used to solve the optimization model including the objective function (1)-(2) and constraints (3)-(9), (14), (19)-(24) in the above steps (1)-(3), to obtain the planned power P_(i) ^(t) of the i^(th) thermal power unit during the dispatching period t and the planned reference power w^(t) _(j) of the j^(th) renewable energy power station during the dispatching period t, in which P_(i) ^(t) is regarded as the planned power of the i^(th) thermal power unit during the dispatching period t and w^(t) _(j) is regarded as the reference power of the j^(th) renewable energy power station during the dispatching period t, to realize a chance-constrained stochastic economic dispatching based on the approximation of the hyperbolic tangent function. 